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I give some scalar field theory calculations on a d-dimensional lune 
of arbitrary angle, evaluating, numerically, the effective action which 
is expressed as a simple quadrature, for conformal coupling. Using 
this, the entanglement and Renyi entropies are computed. Massive 
fields are also considered and a renormalisation to make the (one- 
loop) effective action vanish for infinite mass is suggested and used, 
not entirely successfully. However a universal coefficient is derived 
from the large mass expansion. 

From the deformation of the corresponding lune result, I conjecture 
that the effective action on all manifolds with a simple conical singu- 
larity has an extremum when the singularity disappears. 
For the round sphere, I show how to convert the quadrature form of the 
conformal Laplacian determinant into the more usual sum of Riemann 
C-functions (and log 2). 
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1. Introduction. 



Renyi entropy is a (continuous) extrapolation of Shannon entropy and has 
proved to be useful in various situations. In particular, it has come Tip in the con- 
text of entanglement and a number of explicit calculations have appeared. Casini 
and Huerta, [1], gave some expressions for the coefficient of the (universal) logarith- 
mically divergent term in the entanglement entropy for conformally invariant scalar 
fields on even spheres and, more recently, Klebanov et al, [2], have discussed odd 
spheres. In earlier work, I considered the same situations using, as my workhorse, 
the orbifolded sphere S^/Z^ which introduces a conical singularity of angle 27r/g 
onto the separating submanifold so as to allow the standard contruction of Callan 
and Wilczek, [3] , of the entanglement entropy, 



to be applied. W{q) is the effective action on the orbifolded sphere. 

In the light of the recent activity, I wish to consider Renyi entropies, from the 
same viewpoint but with a different calculational evaluation. Before that, I give 
a further treatment of the computation of operator determinants on spheres and 
extend the analysis to massive fields. Apart from a few underlying manipulations, 
my results are numerical and presented mainly as graphs in sections 5, 6 and 7. 

2. The geometry 

I have referred to the singular manifold as an orbifolded sphere. More exactly, it 
is a fundamental domain for the rotational cyclic action of on S'^. This manifold 
could be termed a periodic lune of angle 'iix /q. I now extend this to a lune of any 
angle, even bigger than 27r. If q is the inverse of an integer this gives a multiple 
covering of the sphere. 

For quantum field theory purposes, I need the spectral properties of the Lapla- 
cian. The conformally invariant eigenvalues on the qf-lune can be taken as the union 
of two sets, Ajv and A^, where, as used before, e.g. [4], 



(1 + qd^) W{q) 



(1) 



Ajv = (a -h -h n2 + . . . + n^f 



\d ^ {.a + qn^+n2 + . . . + n^f 



1 

4' 
1 
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{rii = 0,l,...oo) 



(ni = l,...oo; n2,..,nd = 0, l,...oo). 
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with a= {d—l)/2. These are the eigenvalues on the Neumann and Dirichlet lunes of 
angle 7r/q and follow by standard separation of variables {e.g. Gromes [5], Pockels, 
[6] ) or by algebraic means if q is integral, eg [7] . Any degeneracy is a consequence 
of coincidences. 



3. Even dimensions 



My main interest is in odd dimensions but, for completeness, I repeat a few old 
things. In even dimensions, the effective action generally diverges. Conventionally, 
the coefficient of the associated logarithmic term, say g, provides a universal (i.e. 
regularisation independent) contribution to the entanglement entropy. This coeffi- 
cient, C^i2{q)^ is, essentially, the conformal anomaly and the recipe (1) yields this 
part of the entanglement entropy as the derivative of the contribution to C^i2 of 
the co-dimension 2's worth of conical singularity i.e. , 



9d/2 = -dq{Cd/2{(l) - ^^72(1)) 



(3) 

q=l 

The heat-kernel coefficient, C^i2{q)^ is the value of the relevant ^-function, 
C{s,q), at s = 0. Manipulation reveals that each set of eigenvalues in (2) gives, for 
C{0,q), the (same) sum of two generalised Bernoulli polynomials the properties of 
which then show that the first term in (3) vanishes, i.e. the conformal anomaly on 
the lune has an extremum at the full sphere. This is the essential mathematical 
point of [8] . The conclusion is that the entanglement entropy is just the conformal 
anomaly. 

Turning to the Renyi entropy, the definition is, 

_ nWjl) - Wjl/n) 
^n- , (4) 

where W{q) is the effective action on the periodic g-lune. It can be seen that the 
logarithm coefficient in is essentially just the bracket in (3), evaluated at the 
conventional covering values q = 1/n (although q could be anything). 

It is easy to compute the generalised Bernoulli polynomials by, say, iteration 
and these trivially yield the formulae in Casini and Huerta, [1] , obtained in a differ- 
ent way. For the record, I give the expressions for the Renyi logarithm coefficients, 
g^in) and gio{n), 

{n + l)i79i/' + 79n^ + 23n'^ + 3) 



1814400n7 

(n + l)(1759nS + 1759n^ + 571n^ + lOQn^ + 10) 
239500800n9 
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As a check, setting n to 1 gives, of course, the standard numbers of the classic 
round conformal anomaly.^ 



4. Odd dimensions. Eflfective action. 

For closed, odd-dimensional manifolds, in the absence of zero modes, there is 
no divergence in the effective action using ^-function regularisation and there is no 
logarithm term. In particular, for conformally invariant propagation, the conformal 
anomaly vanishes. In this case, it was suggested early on by Ryu and Takayanagi, 
[11], that the constant term in the entropy (independent of any introduced cut- 
off) could be taken as the universal term. Myers and Sinha, [12], considered this 
viewpoint further in their search for a relevant c-theorem. As pointed out in [13], 
in (^-function regularisation there are no divergences and the entire expression for 
the entropy should be considered universal. The formula (1) still stands, and now 
it is the effective action, that is conformally invariant. 

Before proceeding to novelties, I mention that, in [13], It was shown that the 
effective action, W{q), had an extremum at the singularity-free, round point, q = 1, 
so that the entanglement entropy is just minus the round effective action. Much 
earlier calculations of this quantity could then be drawn upon for analytic forms 
and numerical values. 

Quite recently, Klebanov et al, [2], have evaluated the Renyi entropies for the 
same situation. I will obtain the same results in a more numerical fashion, but one 
that allows extensions to be made. I restrict the discussion to scalars. 

Technically I am assuming, as a working hypothesis, that the effective action 
is given as the mathematically definite quantity. 



where the basic calculational device is the ^-function constructed from the eigen- 
values (2) which I can write as q) = Cni^i + Coi^i To encompass different 
situations, I define a slightly more general object. 



The Bernoulli form of the round conformal anomaly could be avoided by using the integral 
expression derived holographically by Diaz, [9], and by direct manipulation in [10]. 



Wiq)^--C'{0,q), 



(5) 
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where m and u? are d- vectors, so that, 

1 1 

C{s, q) = C(s, a, - I g, 1) + C(s, a + g, - I g, 1) , a={d- l)/2 

^ 1 1 (7) 

= C(s, ojv, - I g, 1) + C(s, a^, 2 ' ^' "^^ ' 

where 1 is a — l)-vector. I refer to the real numbers uj as the parameters^ and 
endeavour to keep them general as long as possible. Also, in the combination (7), 
I require only that + ajy = = 2a;. 

Incidentally, if the field is conformal in [d + 1) dimensions, then a — 0, and 
a, a; u)) reduces to the more elegant Barnes function. In previous evaluations, 
e.g. [4], of the required quantity, C'(0, a, a; u)), an expansion in a sufficed and allowed 
the Barnes function to be employed. Here I wish to work directly with the form 
(6). 

As mentioned in [14], following Candelas and Weinberg, [15], and Minakshisu- 
daram, [16], I can employ the Bessel function form, 

^/ I X a/tt f°° (irexpf— ar) f t \^ ^^"^ ^ , . 

Cs,a,a a; =^/ \ ^ ,A lT-) ^-1/2 «t , 8 

r(s) Jo _exp(-a;iT)) ' 

which exposes the cylinder kernel. 

Adding the N and D parts of the C~function yields the combined expression, 

-1/2 

L/2(«^) 



, , 2^ pc/r cosh ((a;-a)T)/r V-^/\ 
^ ' ^ r(s)7o nti2sinh(a;,T/2)V2«y' ^"'^ 



I{s] 



(9) 



- r(.) ' 

and the object is to evaluate C'(0, a, a | u?) which necessitates a continuation of 
(9) to around s = 0. I adopt the procedure of Candelas and Weinberg, [15], as 
used by Chodos and Myers, [17]. First note that, because the conformal anomaly, 
C(0, a, a I cj), is zero, the integral, /(s), is well behaved."^ Therefore the derivative 
at zero is just /(O). Next remark that the integrand is of the form -7-^*"^"'^ fi'^^) so 
that I{s) can be continued via the complex integral. 



_j_ g7ri(2s — 1 — 



. n^ 9 u /9 ^.-i/2(«T), (10) 

^ J-oo+iA ll^^i 2smh(a;jT/2) V^o;/ 



^ I comment that the value of the C^function at negative integers also vanishes and that the 
derivative at these points is also available numerically as is ^(s) at positive integers. Other 
values would require a complex treatment of Bessel functions. 
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in which s can be set to zero with impunity, so long as d is odd and A hes between 
zero and the first zero of the denominator (which hes on the imaginary axis) i.e. 
A < 27r/maxa;^. (A could be chosen bigger, as long as any poles are allowed for.) 
Doing this yields, 

dr cosh(ct; — a)r coshar 



/(O) = 2 



f dr cosh(cL 



2sinh(a;,T/2) 

1 f°° , ^ cosh (a; — o) T cosh Q!T . 
' dx Re -j , T = X + lA , 



S'' ^ Jo T n?=i sinh(a;,r/2) 

which is a suitable case for numerical treatment for any u), a and a, so long as it 
converges. 

Conformal fields mean a = 1/2. The (i-dimensional g-lune is given by the 
choices, a = {d — l)/2 and u) = (q', 1), so u = {d — 1 + q)/2 and < A < 27r/g. 

I note, for future use, that a massive scalar can be accommodated by setting 
CK^ = 1/4 — A minimal scalar corresponds to a = {d — l)/2 leading to an infra 
red divergence in the integral (11) for large r caused by a zero mode which has now 
to be taken into account. I will not be concerned with this here. 

If q = 1, i.e. the full round sphere, a convenient choice is A = tt, [15], for then 
the real part can be explicitly taken for all d and easy algebra produces, 

(_l)(rf-3)/2 



1(0, d) — Qd-2 / ^(sech a;/2 — sech ~ 

2 Jo X -\- TT 

(_l)(<^-3)/2 



(12) 



2d- 



\j{d) - J{d-2)) 



where 



oo 



J(d)= / dx ^ . (13) 

^ ^ Jo (a;2 + l)cosh'^7ra;/2 ^ ^ 

Equation (12) affords a simpler way of evaluating the effective action on the 
sphere than the sums of Riemann (^-functions resulting from closing the contour in 
(11) in the upper half plane. These sums also follow directly from the manipulation 
of the eigenvalue form of the C^function, the method employed in the earlier calcu- 
lations. A sample number is 7(0,21) = 1.66475568410"^ obtained instantaneously 
from (12).^ As is clear from (12), the values decrease, oscillating with d about zero 
which is, perhaps, not so immediately obvious from the alternative forms. ^ In the 
appendix, I provide a method of deriving the sum of ^-functions form from (12). 



^The Riemann zeta form contains 11 terms. See the appendix. 

different integral expression having the same sign factor was given in [13,18]. An integral 
form, obtained by conformal transformation to a hyperbolic cylinder, is given in [2] for three 
dimensions. 
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For any g, I pursue the simplest path and just compute the integral, 



1(0, d, q) 



2d- 



°° , ^ cothQr/2 coshr/2 
ax Re 



r sinh''"V/2 



(14) 



for any A < 2n/q without further simplification. It is best to plot the results 
as functions of q and the results for the 3d and 5d lunes are given in Fig.l. As 
anticipated, and shown in [13], the effective action has an extremum at g = 1, the 
round sphere. 

Rg.1 . q-lune log det 



o 




-1.5 



-0.5 



0.5 



1.5 



logq 



The proof of this fact in [13] utilised the explicit expression of the effective 
action in terms of Barnes C^functions and some specific properties of generalised 
Bernoulli polynomials. It can be seen here directly from the elementary circum- 
stance that at g = 1 the derivative with respect to q of the integral in (14) is zero. 
There is actually no real part for A = n. 

1 conjecture that this result is a particularly explicit case of a general feature 
of effective actions (functional determinants) on manifolds having a simple conical 
singularity of co-dimension 2 extent such as those, A^^, investigated in [19] and 
[20] , the statement being that the effective action has an extremum when the conical 
singularity fades away (/? 27r). Indeed, it is difficult to imagine otherwise. If the 
lune geometry is deformed from spherical, still maintaining the geodesic embedding 
of the co-dimension 2 submanifold, it is hard to see the extremum disappearing, by 
continuity, and there is nothing to change its location. 
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5. Massive fields 



Because massive fields are not conformally invariant, the significance of the hme 
geometry is not so apparent. However, because the calculation is straightforward, 
1 will investigate an entanglement entropy in the massive case. For the moment, 1 
maintain the assumption (5) for the effective action, which 1 refer to as the 'bare' 
expression. 

Setting o? = 1/4 — /x^, the mass parameter /x measures the deviation from 
conformal coupling and one now has the specialisation of (11), for massive fields on 
the lune, 

^/ , X 1 f°° , cothgr/2 coshar , , 

/(0,d,g,a) = / dxRe . , (15) 

^ Jo T smh'' t/2 

which is, apparently, valid for the entire range of < // < cxo and can again 
be treated purely numerically. The asymptotic limit for large /j,^ can be used as a 
check. Details are given in the next section. 

For the entropy the derivative term, (1), comes into play, the appropriate com- 
bination being, 

1 
2 

_ 1 



e{d, a) = - (7(0, d,l,a) + dg 1(0, d, q, a)) 

dx Re 



9=1 

2 cosh t/2 cosh ar 



T sinh"'T/2 



cosh ar 



sinh''+V/2, 



(16) 



Plots of this entanglement entropy for the 3- and 5-spheres as functions of iJ? 
are shown in Fig. 2. 

An important quantity is the derivative with respect to /U^ at = which is, 
to a sign, the derivative with respect toQ;ata = l/2. 1 find. 



d^&{d, a) 



a=l/2 



1 

¥ 



da; Re 



da; Re 



2coshr/2 
sinh*^"^ t/2 
sinh T — T 



T 



sinh"" t/2. 



=1 r ..... 

2"'io Lsinh''T/2 

/-OO 

= _(_l)(d-l)/2 



cosh'^ a;/2 

m 



7T 



22d-2 r2((d+l)/2) 

For d = 3, this agrees with the value quoted in Klebanov et al ,[21], obtained by an 
unspecified method. The non-zero value shows that this entanglement entropy is 
non-stationary at the UV fixed point, c/ [21]. 
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Rg.2. Massive entanglement entropy 




-0.8 



Q2 0.4 0.6 0.8 1 1.2 
mass squared 



6. Asymptotic behaviour. Renormalisation. 

Information is contained in the behaviour for large mass. The asymptotic ex- 
pansion for the derivative of the (^-function at zero is easily found from the definition 
and is, in its generic form, for odd d, 

C'(0,^')~7^^ E (m^)'/'-'^r(n-rf/2)Q, (17) 

n=0, 1/2,1, ••• 

where the C„ are the short-time expansion heat-kernel coefficients for the operator 
—A -|- y if the actual propagation operator is —A + V + m? . In the present setting, 
it is most convenient algebraically to choose the addition, V , so that —A -|- V 
is conformally invariant in d -\- 1 dimensions. For spheres, and their factors, this 
has the consequence that the heat-kernel expansion terminates, [7]. For spheres, 
V — (((i — l)/2) and the propagation operator is, in our existing notation. 



so that = —m?. 

The coefficients derived in [7] allow the asymptotic form, (17), to be made 
explicit, 

/(0,rf,,,zm)~^ ^~J,kY i?g.((ci-l)/2|g,l)m^ (18) 
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in terms of generalised Bernoulli polynomials. Two examples should suffice 



expressions furnish a useful test of the numerics and the limiting behaviour can be 
seen, roughly, in Fig.l. Because of the termination, the errors in the forms (19) are 
exponentially small, a desirable feature. 

The divergence as m tends to infinity brings into focus the relation between 
the logdet and the effective action since the latter is expected to be zero in the 
infinite mass limit, on physical grounds. Hence I could take the right-hand side 
of (18) over to the left hand side and consider this to be a finitely renormalised 
(minus) logdet. This is equivalent to the renormalisation procedure advocated by 
De Witt, [22], and amounts to subtracting from the propagation heat-kernel, K{t), 
sufficient terms of its short-time expansion so as to render the formal expression 
for the effective action. 



finite. (See [22], equn. (14.75).) 

The problem with this is that the so renormalised effective action has to be 
evaluated at the conformal point, = 1/4 i.e. at = —1/4 where it is non-zero 
and, worse, formally imaginary. Hence I will make a further, ad hoc hypothesis 
and require the subtraction term to vanish at this point, in addition to having the 
asymptotic form, (18). This can be accomplished by shifting the mass to = 
1/4 — which vanishes at the conformal value. The term to be subtracted from 
C'{0, q) is then in three dimensions. 



which for the three dimensional qf-lune is to be compared with (19) with Cg = 27r^/g 
and = 27r^(g^ — l)/3g. The renormalised massive effective action for this case is 
thus defined as. 







(20) 
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where /(O) is given by (15) with ct^ = 1/4 — /i^. A few plots of this effective action 
are presented in Fig. 3. 



Rg.3. renormalised effective action 




ns 
CP 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
mass squared 

This renormahsation procedure can at best be considered makeshift. An ad- 
vantage is that the vanishing at the conformal point means that aU the preceding 
conformal arguments hold true, including the extremum at q = 1. However, the 
infinite n behaviour is not exponentially vanishing. Because of this I do not dis- 
cuss the higher spheres, but I will give the resulting form of the entropy for three 
dimensions. The renormahsation modifies the bare entanglement entropy &{d, a). 



®ren(3,a) = ©(3, tt) + ^ 



(a 



1/4-^' 



for d — 3 and a plot, against /u, is given in Fig. 4. 

Rg.4. 'Massive entanglement entropy 
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7. Renyi entropies 



In Fig. 5, I plot ttie conformal Renyi entropy, S^, (4), versus a continuous 
n = 1/q for tlie 3- and 5-spheres. 

Rg.5. Renyi entropy 



a. 



m -0.5 



T 1 1 1 1 1 1 r- 
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J I I I I I I L_ 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

n 

Fig. 6 is for the 3-sphere only, and shows the variation of the renormahsed 
Renyi entropy with mass /i. Using (20), this is given, for d = 3, hj 



(21) 



Fig. 6. massive Renyi enlropy 




mass 

The coefficient of the linear asymptotic behaviour of with mass (cancelled 
in (21)) is expected to be a universal constant (see [23]). For arbitrary dimensions, 
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the form, (18), aUows one to find this constant as the coefficient of the highest power 
of m (or /i), which is m^~'^. Easy algebra yields for this coefficient, 

"^"^"^^ ^ 12(d-2)! n ' 

for the present toy model. This should be compared with the similar Gaussian 
expression obtained in [23] , cf also Calabrese and Cardy, [24] . 
Setting n = 1 gives the entanglement entropy constant, 

r,(<i) = (-l)(-W2_^, (22) 

which could be compared with the formula of Hertzberg and Wilczek, [25] , derived 
in basically the same way, but for a different geometry. Their expression, on taking 
the dividing surface to be a (d — 2)-sphere, turns out to be identical to (22) showing 
that the curvature makes no difference to the result which is, not surprisingly, due 
just to the local conical singularity. 

Hertzberg, [26] , has analysed the dependence of the entanglement entropy on 
the mass for an interacting field theory. To leading order, the result is the replace- 
ment of the bare mass in the free field expression by its renormalised value. 

The (finite number of) lower powers of the mass in the asymptotic form are 
curvature effects. This is in concordance with the observations of Lewkowycz, Myers 
and Smolkin, [27], who generalise the geometry of [25] to a waveguide of curved cross 
section. 



8. Discussion 



I have presented a more numerical approach to the computation of the confor- 
mal effective action, and related quantities, on an odd-dimensional lune of arbitrary 
angle. For the full, round sphere this leads to yet another expression for the Lapla- 
cian determinant. In particular I thence obtain the older form of a sum of Riemann 
^-functions at odd arguments. 

The extension to a massive field is made and I attempt to renormalise the effec- 
tive action so that it vanishes for infinite mass. The renormalisation implemented 
in section 6, which forces the effective action to vanish also at the conformal point, 
is suspect for small values of the mass, however I could not devise anything better. 

Going beyond the lune, the numerical quadrature expression, (11), allows the 
log determinants on other, special spherical regions to be computed by an appro- 
priate choice of the parameters, uj. For example, the conformal log determinant 
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on the periodic tesseract (Schlafli symbol [3,3,4] with u? = (4, 6, 8)) is 0.579184 and 
will be discussed at another time. Related is the computation of Barnes' multiple 
Gamma function, depending on a and w, which is quite hard to find. 



Appendix. Expression in terms of Riemann function. 

For those who like such things, I outline a way of transforming the integral 
expression (12) into the alternative form of the effective action on the full sphere 
as a sum of Riemann ^-functions evaluated at positive odd integers. I also include 
some ancillary mathematical remarks. 

The basic notion is to expand the powers of sech x/2 in multiple even derivatives 
of secha;/2. For odd d = 2r + 1 this is clearly possible ^ and I define coefficients, 
by, 

sech^'"+'a;/2 = J2^p ^ secha;/2, (23) 

p=0 

so that, from (13), one requires the integral. 

Kip) = / dx — — —r- sech7ra;/2 

1 [°° d'^P 1 

7^^::2T sech TTx/2 (24) 



TT^P Jq dx'^P (1 + a;2) 
(2p)! r sin(2p + 1)6 

TT^P 



f°° , sin(2p + l)6' , 

io (1 + X2)2P+1 '^^^^^/^' 



where I have used Liouville's form of the multiple derivative (see Gregory, [29], 
p. 17, no. (20)) and 9 — tan~^ 1/x. The oddness of the dimension (evenness of the 
derivative) means that sin(2p + 1)6' = (-1)^ cos ((2p + 1) tan"^ x) and at this point 
I recall an integral form for the Riemann (^-function due to Jensen which, at odd 
arguments, yields, ^ 

2'^P r°° 1 cos ((2p + l)tan-ia;) 
Ci?(2p + 1) = ;l _ 2-2p (l + a;2)2p+i cosh7ra;/2 " 



^Something related, actually the inverse, is performed by Stern, [28]. 

^ This formula is obtained in Lindelof , [30] , using a variant of Plana- Abel summation applied to 
the function and it may be a more uniform complex treatment of the present identities can 
be found. 
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It is therefore seen that the quantity J{d), from which the derivative at zero, 
1(0, d), can be found using (12), is given by. 



J{2r + l) = J2^;K{p) 

p=0 

= E(-l)'' Ch(2p + 1) + log 2 , 

p=l 

which exhibits the usual structure.^ The end value J(l), (associated with the pole 
of the Riemann ^-function) is a special case and {e.g. Gregory, [29], p. 496 example 



Jo (a^ + l)cosh7ra;/2 
Incidentally, a more general integral is, [31], Table 97, No.4, 



(x^ + a^) cosh7rx/2 \ 4 

1 ^(a-l)/2 



2 



+2; 



in terms of Stirling's /3 function (e.g'. Nielsen, [32], pp. 16, 181) and Legendre's inte- 
gral form. 

There is no simple expression for the rational numbers Sp but they can be 
rapidly computed from the recursion, somewhat similar to that for generalised Euler 
numbers, 

- ~ r(2r - 1) ^f-^ ' 

together with SL\ = ^ = ^r+i the initial value £q = 1. 

Expressions for the effective action in this particular form can be found listed 
in [33] up to d = 11. Out of interest I present the vector of the coefficients of 
Ciii2p+1)/2VP for p from 1 to 6 (i.e. d = 13), 

2385868 21914 6932 7174 ' 

, , , , -4774, -4095 

17325 ' 135 ' 15 ' 3 ' 

The coefficient of log 2 is 21. 2"^^ 

^ Instead of the ^-functions at odd positive numbers, the derivatives at negative even integers 
can be used. This results in a certain cosmetic improvement. 
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The numerical evaluation of these expressions can easily be automated, but is 
not so efficient as the earlier one using quadratures. 
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